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ABSTRACT 


Under typical conditions in the subtropical thermocline 


double-diffusion occurs in the form of salt fingering 


(warm, salty fluid overlies cold, fresh fluid). The 


formation of staircases in the thermohaline structure of 


the ocean has been observed since the late 1960s, with 


recent field data collected within the tropical Atlantic 


displaying staircases with high-gradient interfaces 
characterized by a unique spatial orientation determined by 
background temperature and salinity. Competing theories 
have been proposed to explain the dynamics of these 


staircases; however, the origin of the staircases and the 


mechanism that controls final equilibrium remains poorly 


understood. This thesis examines staircase development in 


the tropical Atlantic. Incorporating double diffusion using 


the flux-gradient formulation of Radko and Smith, 


staircases are numerically simulated to resolve the 
controversial aspects of the  staircas theories. The 
staircase simulations are critically evaluated against the 


C-SALT experiment observations, and the conditions for 


their formation and explanation of the fully equilibrated 


state are resolved and explained. The effect of acoustic 


propagation through the three-dimensional modeled 


staircases is evaluated to determine the impact of these 


large staircase areas on various frequencies and depths. 
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LD. INTRODUCTION 


A. DOUBLE DIFFUSION 
1. History 


Double diffusion is a phenomenon that occurs due to 


the differing molecular diffusivities of two properties 
(e.g.; temperature and salinity), within a stably 
stratified fluid. Double-diffusion operates in 


configurations in which these properties have an opposing 


effect on the stability of the water column. 


Double diffusion is a relatively recent discovery in 


the realm of science (Stern 1960), although W. Stanley 


Jevons, an English economist, came very close to 


discovering it in 1857 when he created sugar fingers during 


an experiment. Although Jevons understood that heat 


diffused faster than the dissolved sugar, his statement 


that “all convection should occur as minute streamlets” 


failed to capture the physics behind th phenomenon 


tschma.tt 1995): 


Whilst others would go on to attempt to explain the 
phenomena, including the English physicist Lord Rayleigh, 
the subject was correctly explained only in the mid-20th 
century by Melvin Stern (1960). Stern’s discovery was 
motivated by an experiment by his colleagues, Stommel, 
Arons and Blanchard (1956), who developed the Perpetual 


Salt Fountain. This was an experiment initially designed to 


test pressure at the base of an ocean, by inserting a 
narrow heat-conducting pipe vertically to draw water to the 


surface. Due to the lower salinity of the basal water, once 


at thermal equilibrium it would continue to rise towards 
i, 


the surface due to being lighter than its surrounds. While 
investigating this idea, the concept that placing a spigot 
on the top of the tube would cause the water to run 


continuously was introduced, thus the “Perpetual Salt 


Fountain” was born. Stern took this experiment to the next 


step, realizing that salt fountains occur naturally in the 


ocean, as a result of heat conducting 100 times faster than 
the diffusion of salt in the = ocean. Hence thermal 
equilibrium of a parcel of water with its surroundings can 
occur faster than chemical equilibrium, resulting in an 


ocean containing parcels with differing salinities adjacent 


to each other devoid of a salt fountain pipe. Different 
salinities result in different density gradients, and thus 


motion dubbed as double diffusive convection. 


2. Salt Fingers 


There are two types of double diffusive convection in 


the ocean. The first is when cold, fresh water overlies 


salty, warm water. If a warm salty parcel is displaced 


upward, due to the difference in molecular diffusivity it 


cools much quicker then it loses salinity, thus becoming 
denser than its surroundings and sinking back down. As it 


Sinks it regains heat quickly, but is now less saline than 


its surroundings, thus rising again, this time further than 


it did previously. This process can continue to occur, 


slowly mixing the water column. This is known as diffusive 


convection. The other type is salt fingering, which occurs 


when warm, salty water overlies cold, fresh water. If the 


warm, salty parcel is displaced downward, as in the 


previous example it will lose heat quicker than salinity, 


becoming denser than its surroundings, and thus continuing 


2 


to sink. These parcels appear as thin elongated streams, 


only 1 to 5 centimeters wide, and tens of centimeters long. 


Salt Fingers 


Warm, 
Salty 


pat = 
xchange 


Cool, 
Fresh by diffusic 


Sat carried 


downin —" 
fingers 
Profiles S(z) T{z) P(Z) 
Final Initial Final} !nitial ‘ 
» Final) ;!nitial 
i 
! i | | 
Salt Flux Heat Flux Density Flux 
Figure 1. Schematic diagram illustrating the physical 


mechanisms of salt fingering (From Ruddick and 
Kerr, 2003) 


32 Thermohaline Staircases 


In the ocean, vertical stratification is often 


characterized by the appearance of well-mixed layers 
separated by fine interfaces with large temperature and 
Salinity gradients, known as thermohaline staircases. The 
interfaces of these staircases are occupied by salt 


fingers, which cause a net unstable buoyancy flux that 


stirs the layers (Stern and Turner 1969). This stirring of 


the layers can result in vertical growth of the layer 


through destruction of the interface and merging of the 


layers. 


There is a large abundance of these staircase regions 


throughout the world, and the increased temperature and 
salinity fluxes in these staircases (Schmitt 2005) suggest 


that they play an important role in global ocean mixing. 


4. Theoretical Background 


The density ratio controls the intensity of salt 
fingering and is defined by: 

_oAT 

P BAS (1) 

where, @ is the thermal expansion coefficient and B is the 

saline expansion coefficient. If the density ratio is less 

than one, then the density gradient is unstable and will 


overturn. High values of density ratio correspond to weak 


double diffusive mixing, whereas density ratios that are 


slightly above unity are associated with vigorous 


fingering. 


Another important parameter in the formation of salt 


fingers is the flux ratio, 


where F, and F, are the molecular fluxes of heat and salt, 


respectively. This value must theoretically be below one to 


satisfy the net transport of density “up gradient.” 


B. C-SALT 


In 1985, a field program by the name of Caribbean- 
Sheets and Layer Transects (C-SALT) was carried out in the 


tropical North Atlantic Ocean east of Barbados to survey 


the thermohaline staircases in the area, developed due to 


vertical mixing by salt fingers. The survey area lies 
between 150m and 750m depth, at the confluence of the 
highly saline Subtropical Underwater and the fresher 


Antarctic Intermediate water, and this gives rise to a very 


strong destabilizing salinity gradient which can drive salt 


fingers (Schmitt 1987). The temperature profiles in the 


survey were classified as layered, weakly layered or 


nonlayered, with an area within 10° and 14°N, and 52° and 


57°W having the thickest mixed layers, largest number of 
layers, and largest temperature differences, as seen in 


Figure 2. 
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Figure 2. Area of C-SALT survey (From Schmitt, 1987) 


It is this area that becomes the focus herein as we 


attempt to analyze, simulate and theorize the processes 
that occur within this staircase region, and conduct 


acoustic modeling through the staircases. 


II. LATERAL COHERENCE, MERGING, AND ALIGNMENT OF 
THERMOHALINE STAIRCASES 


A. LATERAL COHERENCE 
1. C-SALT Observations 


An interesting feature of the C-Salt observations was 
lateral coherence of the staircase layers over many 
kilometers by examining profiles at a variety of 
separations. The shortest separation between consecutive 
profiles was 1.3 km using a CTD tow-yo, covering a total 


distance of 33 km. The results of this are shown in Figure 


3, whereby individual layers can be traced over the entire 


hydrographic section. 


TEMPERATURE 


Bs 
a 
— 
TEMPERATURE 
Figure 3. Three-dimensional view of CTD tow-yo (rear 


profile to the south) as the ship steamed to the 
north. 1.3 km separation between profiles, 33 km 
total distance covered (From Schmitt et al. 1987). 


This lateral coherence was not limited to 33 km, as it 


was also measured in individual layers over a distance of 


more than 100 km using aerial dropped bathythermographs 
(XBT) spaced every 5.5 km in an east to west pattern 


(Figure 4). 
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Figure 4. XBT temperature profiles from east to west, with 
5.5 km separation between drops and total distance 
of 220 km. The profile at far left is correct, each 

subsequent profile is offset by 1.6 °c (From 
Schmitt 1987). 


Over longer distances there were Significant 
temperature and salinity changes within the layers, however 
an examination of the temperature and salinity 
characteristics of these layers produced a _ surprising 
pattern where temperature and salinity values are grouped 


in distinct lines, suggesting layers get warmer, saltier 


and denser to the north and west (Schmitt 1988). Figure 5 


displays this unique grouping, which demonstrates that the 
well-defined slope of the layers crosses isopycnals. The 
horizontal density ratio, when averaged over the ten 


principle layers, was 0.85 (1+0.02) in the spring, and 0.84 


(40.03) in the fall. Among other things, this value aided 
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in the realization that salt fingers were the dominant 


process in maintaining C-SALT staircases, as the value of 


the lateral density ratio was distinctly different from 


vertical mixing by small-scale turbulence (R,=16), and 


isopycnal eddy-stirring (1.0) (Schmitt 1988). 
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Figure 5. Potential temperature—-salinity values of layers 


from CTD stations during C-SALT. Solid circles are 
from layers greater than 10 m thick, open circles 
are from layers 5-10 m thick. Hashed lines are 
temperature-salinity relationship obtained from 
stations at northwest and southeast corners of the 
survey. Solid lines represent isopycnal surfaces 
(From Schmitt 1987). 
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B. MERGING 
d. Theory 


It is difficult to gauge the evolution of layers in 


the C-SALT observations, short of saying that some 
staircases have persisted over 15 years, having again been 
observed in 2001 (Schmitt 2005). The evolution and long 


term persistence of these staircases can only be 


speculated, and the question remains as to whether they 
form over decades or centuries, and what, if anything, 
causes a collapse in the staircase. The theory behind 
staircase layer merging is explored here, and validated 


using a numerical model discussed later. 


It has been proposed, and confirmed by direct 


numerical simulation (DNS), that layer formation occurs as 
a result of the instability of the heat/salt flux ratio, Y, 
that occurs when the flux ratio is a decreasing function of 


density ratio (Radko 2003). Further to this proposal is the 


theory that relatively thin layers within the staircase 


merge until the thickness of the layers exceeds a critical 


value, H with equilibrium occurring as a result of 


o7 
inhomogeneity of the convecting layers stabilizing the 


staircase (Radko 2005). 


In Radko (2007), a merging theorem was developed which 


states: 


If the buoyancy flux decreases with the buoyancy 
variation across the step, then the staircase is affected 
by the B-instability, and the spontaneous merger of layers 
occurs through gradual erosion of weak interfaces. If, 


however, fluxes increase with step height, then coarsening 
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of a staircase may result from drift and collision of the 


adjacent interfaces. [Radko, 2007]. 


In Chapter III, we analyze the merging events in the 
numerical model, using the initial temperature of a layer 


at commencement of a merging event and the temperature 


variation in the layer during the event. In all simulations 


analyzed in this study, mergers occur following the B- 
merger scenario. 
Cc. ALIGNMENT WITH BACKGROUND FLOW 


Another area of research extending from the C-SALT 


experiment that was explored further by Wall (2007) was the 


alignment of the slope of staircase interfaces with the 
background large scale isotherms within the staircase 
region. Wall first began by considering a mixed layer 
sandwiched between two interfaces as shown in Figure 6, and 


starting with the conservation law for temperature and 


salinity in a mixed layer (F,, and Ff, are the temperature 


flux at the top and bottom of the layer, respectively). 


1. 


ia 


T1 


Figure 6. Basic two layer interface (from Wall 2007) 


By adopting a co-ordinate system that is in the same 


direction as the flow, and using the Radko (2005) 


hypothesis that as a staircase evolves in time to a quasi- 


equilibrium state, the interfacial flux ratio approaches 


its minimum value, Yi; Wall arrived .at. the following: 


ae / 03) = (En. Zz Fad = Vo (3) 


(Si9./ 0%); Eg =i) 


Wall rearranged the left hand side of this equation to be 
(—), and determined this value from the temperature- 


salinity plot similar to that shown in Figure 5 for 
stations in the C-SALT area exhibiting strong staircases. 
The temperature-salinity variation along the mixed layer 
was derived in terms of Cartesian derivatives and after a 


series of derivations Wall produced: 
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where @, and Q@, are the slope of the background temperature 


and salinity, respectively, a@ is the slope of the layer 
interface, and assuming R,>Y,,, yields: 
tea ~ tgdl, (6) 


Wall thus hypothesized that the interface should align with 


isotherms but not necessarily isohalines. An observational 


test was carried out on the C-SALT staircases discussed 
above, uSing a numerical model based on the GFDL MOM3 code 


(Pacanowski and Griffies 1999) to determine th direction 


of the mean background flow. The results of this 
observational test showed a weak relationship in the 
direction of background flow between the isotherm and the 


slope of the layer (much less so in the case of 


isohalines). In the direction normal to background flow, 


the relationship between the isotherm and slope of the 


layer disappeared. This idea is further explored in the 
following chapter, in order to attain a stronger 


confirmation of this feature shown by Wall. 


1.3 
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IIIT. SIMULATION OF STAIRCASES—-EVOLUTION AND 
ALIGNMENT 


Ay MODEL SETUP 


5 Massachusetts Institute of Technology General 
Circulation Model (MITgcm) 


In order to simulate the staircases seen during the C- 
SALT field program, the MITgcm numerical model was chosen. 


The desirable features of this coupled model include, but 


are not limited to, the non-hydrostatic capability that 


allows the model to simulate fine, meso, and synoptic scale 


processes, and the ability to perform on a wide variety of 


computational platforms (Adcroft et al. 2012). The 


flexibility in building the model with varying parameters 


was ideal for the purpose of this investigation, where many 


parameters such as the horizontal and vertical diffusivity 


rates are not known, but can be slowly adjusted to produce 


simulations that reflect conditions in nature. 


Model runs were performed on the Department of Defense 


High Performance Computing Modernization Program, Cray XE6 


Open Research System (Chugach), located at ERDC, Vicksburg, 


Mississippi. 128 processors were used in coarse resolution 


grid simulations encompassing the entire ocean basin. 


During early model initialization and simulation a 


coarse uniform 1 degree grid resolution was used over th 


entire ocean basin. Once this model had reached a period of 


perceived equilibrium at time = 18 yrs, a high resolution 


(Ax = Ay = 10 km) interpolation was carried out in the area 


of staircases, with the remainder of the ocean basin 


staying at coarse grid resolution. The simulation used Az=l 
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m and a Cartesian grid. Initial runs involved the use of 


the Gent McWilliams/Redi SGS eddy parameterization scheme 
(GMredi) which mixes tracer properties along isentropes 


(Redi 1982), and adiabatically rearranges tracers through 


an advective flux where the advecting flow is a function of 


the slope of the isentropic surface (Gent and McWilliams 
1990). The GMRedi scheme was implemented using the tapering 


scheme of Large et al. (1997); however, it was found that 


over time the vertical temperature salinity profiles were 


smoothed to such an extent that the staircases were 


removed, an unexplained result in the simulation in direct 


opposition to what we know happens in the ocean. Thus the 


GMredi scheme was not used in most simulations. 


Trial simulations were also carried out using various 


values of biharmonic diffusion coefficients, in order to 


accurately account for both mixing as a result of 


turbulence and double diffusion. The biharmonic diffusion 


coefficient was ultimately parameterized at 1*10° m‘/s over 


all 1000 vertical layers. 


Various parameter adjustments were made throughout the 


course of model xperimentation, and these are discussed 


later in the chapter. 


2. Initial Datasets 


At model initialization, temperature and salinity 


datasets were generated to provide a density ratio pattern 
over the tropical north Atlantic qualitatively resembling 


observations during C-SALT. The initial modeling of 
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temperature and salinity used 


the surface and 1000m is 


values in Figure 7. 
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Temperature and Salinity at the surface and at 
and Density Ratio, at initialization. 


1000m depth, 
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3: Model runs performed without Double Diffusion and 
with Constant Flux Ratio 


In order to gauge the effect of double diffusion on 


the simulation, several model runs were conducted using 


different K profile parameterization (KPP) setups. The KPP 


scheme within MITgcom is a product of Large et al. (1994), 


and it unifies the treatment of a variety of unresolved 
processes in vertical mixing (Adcroft et al. 2012). The 
functionality to switch double diffusion on and off within 
the KPP scheme allows the ability to display a profile 
within the expected staircase region with and without 


double diffusive mixing, so as to se th effect double 


diffusion (and in this case, salt fingers) has on the 


development of staircases. Another useful feature of the 
KPP scheme in MITgcm is the ability to adjust the flux 
ratio. The flux ratio model used for all modeling is 


suggested by Radko and Smith (2012), who used high 


resolution DNS data to obtain the following algebraic 


expression: 


Y = a, exp (-B,B.) + C, (7) 


where (ay, b,, C,) =<(2.5709, 2.541:3, 02-5128). 2 To illustrate the 


effects of using a variable flux ratio, another simulation 


was also conducted with a constant flux ratio of 0.5128. It 
is suggested in Radko and Smith that for the step-like 


patterns of thermohaline staircases to occur, the flux 


ratio must decrease with an increasing density ratio, in 


the oceanographically relevant range (1<R, <3): We put this 


theory to the test by comparing these two simulations with 
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the expectation that with a constant flux profile a 
smoothed vertical temperature profile would be evident, as 


opposed to a flux ratio utilizing Equation 10, where a 


stepped vertical temperature profile would be expected. 
Figures 8 through 10 show a comparison of the varying KPP 
examples discussed here at time, t#18 years, and it can be 
seen that the ideal KPP set up for growth and evolution of 


staircases similar to those observed in C-SALT is using the 


flux ratio at Equation 7 with double diffusion as a 


dominant vertical mixing process. Profiles are taken at the 


1 1 
location rae and i 


representative of the stratification in the entire domain. 


P but the patterns are 
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Figure 8. Vertical Temperature profile with constant flux 
ratio, y=0.5128. 
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Figure 9. Vertical Temperature profile with the effects of 
double diffusion removed from the model. 
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Figure 10. Vertical Temperature profile with flux ratio from 
Equation 10 and active double diffusion. 
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B. MODEL RESULTS 
1. Effects of Vertical Eddy Diffusivities 


Vertical eddy diffusivities due to small-scale 


turbulence, K, and K,, could be specified within MITgcm 


and typical oceanic values are of the order 10° m/s. In 
order to examine the impact of turbulent diffusivities on 


thermohaline layering, the model was assigned a series of 


K, and K, values ranging from 0 to 10* m’/s , and then run 


to t#18 years on each simulation. A Matlab code was 


developed to compare the staircase distributions with a two 


dimensional map of the density ratio, Roy between 750 and 


850 m. To identify staircase distributions, the average 
density ratio over the 100 m vertical section was 
calculated, and areas where density ratio changes within 


this area wer greater than 10% of th averag wer 


identified: 


max R. —- minR 
=) eS AR (8) 
: 1000 


where R, is 10% of the average density ratio and AR, is the 


change in density ratio. If the individual density ratio 
change between 750 and 850 m was greater than 10% of the 


average density ratio change, and more than a third of the 


vertical profile (between 750 and 850 m) at individual 


points exhibited this characteristic, than staircases were 


deemed to exist at this location. 


The results are displayed in Figures 11 through 14, 
where the white contours indicate the presence of 


thermohaline staircases. Alongside the staircase area is a 


graphical representation of the flux ratio dependence on 
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the density ratio at the specified eddy diffusivity. As per 


the suggestion by Radko and Smith (2012), thermohaline 


staircases are prevalent when the flux ratio decreases, but 


are almost non-existent when the flux ratio increases with 


increasing density ratio. It is evident that when K,=2.5* 
10° m*/s and K,=2.5*10° m/s, the numerical model closely 
replicates the observed density ratios and staircase 
formation areas of the C-SALT experiment, and at these 
regions the flux ratio is also decreasing. Other parameters 
that were adjusted to best replicate C-SALT observations 


are in Table l. 


Parameter Value 
Vertical Eddy Viscosity 1*10* m/s 
Lateral Eddy Viscosity 2*10° m/s 
Lateral Eddy Diffusivity 1*107 m/s 
Table 1. Parameter values in the numerical model. 
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Figure 11. Total flux ratio versus density ratio (left) and 


density ratio (right) for K,=T,=0 m/s. Upon 
analysis staircases are evident throughout the 
region. 


Longitude (deg) 
Figure 12. Total flux ratio versus density ratio (left) and 
density ratio (right) tor K, =K,=10° m/s. 


Staircases are largely prevalent throughout the 
region as indicated by contours. 
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Figure 13. Total flux ratio versus density ratio (left) and 
density ratio (right) for Ku=K,=25*10" m/s. 
Staircases are evident a region similar to that 
observed in the C-SALT experiment. 
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Figure 14. Total flux ratio versus density ratio (left) and 
density ratio (right) for K.=K.210" m/s. 
Staircases are almost completely absent at this 
vertical eddy diffusivity. 
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2. Mechanism of Merging and Growth Rates 


The model output was examined from initiation to 


perceived equilibrium. The images in Figure 15 through 17 


represent a three dimensional visualization of the 


Staircase looking north east from 50 West and 8 North 


between depths of 500 to 800 m shown at initialization of 


the model, an early point in the merging stage of the 
staircase, and fully matured merged layers. This three 
dimensional rendering is produced by VAPOR 


(www.vapor.ucar.edu), a product of the Computational 


Information Systems Laboratory at the National Center for 


Atmospheric Research. 


Figure 15. Three-dimensional visualization at t=0 secs. 
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Figure 16. Three-dimensional visualization at t=10 years. 


Staircases are evident at 
commenced growth 


this time and layers have 


through merging. 


Figure 17. Three-dimensional visualization at t=18.6 years 


with evident staircases merged to equilibrium. 
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The evolution of the temperature profile in time is 


shown in Figures 18 and 19, from the initialization (t=0), 


the time t=5.4 years at which flux ratio instability 


modifies the initial profile, the time t=9.2 years where 


numerous density inversions are evidenced within the 
profile, to the point where a well-defined staircase is 
evident at t=10.4 years, and then several time steps 
illustrating the merging of the steps through to perceived 


equilibrium at t=18.6 years. Profiles are taken at the 


or 50° W, 8 WN on a global 


1 1 
location are and are 


grid. 
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Figure 18. Formation and evolution of layers in the 
numerical model shown at (a-d) t=0, 5.4 , 9.2, 10.4 
years. 
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Evolution of layers in the numerical model shown 


18.6 years. 


revealing the merging events over 


the 


interface growth at the expense of the 


the B-merger as discussed in Chapter II). 


of time 
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Figure 20. Space-time diagram of the numerically modeled 
temperature gradient (c/m). Layers merge when the 


relatively strong interfaces grow further at the 
expense of weaker interfaces, which decay and 
eventually disappear (B-merger scenario). 


Essential insights into the layer-merging process 
brought by theoretical studies (Huppert 1971 and Radko 
2005), is that the merging effect is a form of linear 
instability. Thus perturbations for a quasi-steady state 
are expected to grow exponentially. Our simulations offer a 


convenient opportunity to test these theories. Figures 21 


and 22 show two separate merging events selected from the 
DNS over the time frame of the merge. The growth rate is 
exponential for both examples, which supports the theory 


merging as a result of linear instabilities. 
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7s Lateral Coherence of Modeled Layers 


Another featur of C-SALT that was also evidenced in 


the MITgcm simulation was lateral coherence of the layers 


over many hundreds of kilometers. The lateral extent of 


these layers within the simulation measured greater than 


300 km with minor spatial changes in temperature and 


salinity within the layers as was also evidenced in C-SALT. 


Temperature and salinity data from three individual layers 


is plotted in Figure 23, and the result clearly portrays 


the layer temperature and salinity grouping aligned in a 


straight line with the layer slope crossing isopycnals, as 


it did during C-SALT (Figure 5). It is noted that this 
pattern became more pronounced as horizontal eddy 


diffusivities, K,T and K,S decreased by 0(10), from 1000 to 


100 m/s. Th average density ratio was defined by the 


temperature and salinity changes and was calculated to be 

2.01, and the resultant theoretical flux ratio using the 
equation of Radko and Smith is 0.53. The line of this slope 
is shown in Figure 23 below each set of layer points to 
show the correlation between the flux ratio and lateral 


density ratio. 
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Figure 23. Temperature-salinity values taken from three 


individual layers in the model. Each color 
represents a separate layer. 


The lateral coherence is clearly displayed in Figures 


24 and 25, which show a series of vertical temperature 


profiles from the numerical model taken in the C-SALT area 


at horizontal separation of 10 km in the east-west 
direction, covering a total distance of 290 km. The 
temperature scale is correct for the far left profile, with 
each consecutive profile offset by 1.6 °C. Individual 
layers can be traced along the entire distance, displaying 
a very clear lateral coherence, which extends over 300 km. 
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Figure 24. Series of offset temperature profiles from the DNS 

in an east to west section with 10 km separation. The 
right-most profile lies to the east of the others; 

the total distance covered is 290km. The temperature 

scale is correct for the left-most profile and each 


subsequent profile is offset 1.6 C.. 
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Figure 25. Temperature profiles from Figure 24, focused at 
the depth of staircase formation. 
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4. Alignment of Interface Slopes with 
Flow 


To attempt to correlate the slope of th 


Background 


interfac 


with the background flow as discussed in Chapter II, the 


background flow must first be found. A Matlab code was run 


to identify the lateral temperature and salinity, 


along the direction of the flow using the 


formulation, where u and v are velocity in 


respectively, in a Cartesian grid: 


jay ceeoiaeey yen 
vy 
Uu V 
8-38, +5, 4 
V| v 


T, and S,, 
following 


x and wy, 


(10) 


The slope of the isotherm and isohalines can then be 


defined as: 


qT 
tg(a.,) = 2 
ek 
S) 
tg(O,) = Ss 


Zz 


(11) 


The mean background flow was found to be predominantly 


West South West in the DNS throughout the C-SALT area of 


interest. The intent was to determin th slope of the 


interface by identifying the layer, and then 


determining 


the gradient change. Unfortunately, as can be viewed in the 


lateral coherence figures, the depth change of 
was exceptionally small over several hundred 
(layers are almost completely horizontal). This 


a zero degree angle for the interface, th 


the layers 
kilometers 
results in 


refor the 


following explanation of the continued analysis is provided 


for future research purposes. 
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Using Wall’s formulation for the interface slope 


(Equation 7), the theoretical value can be compared with 


the numerical simulation. Taking Eguations 9 and 10 and 


shifting 90 degrees to the right to evaluate normal flow 


gives: 
T= ae SY 
yo iy 
(12) 
Ss. = pedis 
| | 


where the numerator in Equations 13 and 14 is now given to 


be T, and S,, respectively. It is proposed that this can 


also be compared with the interface angle to demonstrate a 
lack of correlation between interface slope and isotherms 


or 1isohalines normal to the direction of flow. 
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IV. STAIRCASE EFFECTS ON ACOUSTIC PROPOGATION 


A. BACKGROUND 


The effect of thermohaline staircases has previously 


been studied, both within the C-SALT area and in other 


areas of the ocean. Chin-Bing et al. (1994) used data that 
was collected from the northeast coast of South America 
adjacent to the C-SALT area in 1985. These staircases were 
formed as a result of salt-fingering. The study focused on 
the effect of the sound velocity profile as a result of the 


staircases, and resultant acoustic propagation simulations. 


The study only investigated the effect at 50 Hz, which 


creates an acoustic wavelength of similar size to the step 


size. A conclusion from the study was that for significant 
differences in prediction to be observed both the source 


and receiver needs to be placed within the depth range of 


the staircases. 


Wall (2007) built on this investigation of acoustic 
propagation using a parabolic equation (PE) model for a 
stepped profile taken from C-SALT. The PE model was run 
using two different source depths (one within and one below 
the staircase region), and three different source 
frequencies. Contrary to Chin-Bing’s findings, Wall found 


that the source and receiver did not need to be at the same 


depth as the stepped layer for significant differences to 


be observed between a smooth and stepped profile acoustic 


prediction. This finding was similar to one found in the 


past by Wilson (2007), who studied propagation through a 


stepped structure resultant from diffusive convection. Wall 


also found that as the frequency of the source is increased 


37 


the effect of the staircase on the acoustic prediction is 
amplified. A limitation to these observations was 


frequency, where 1200Hz was the highest frequency studied. 


This study continues to build on Wall’s observation, 


using a different model (Bellhop) to compare past results, 


and increasing the frequency to determine if the decreas 

in wavelength continues to have an impact on difference in 
acoustic prediction between smooth and stepped profiles 
taken from C-SALT. The analysis of the source within the 
staircase is investigated further than the previously 
mentioned studies, such that the source is set at a depth 


corresponding to a prominent interface, and also at the 


center of a prominent mixed layer to determine whether the 


positioning of the source within the staircase region plays 


a pivotal role in the propagation. Based on the findings of 


lateral staircase alignment in earlier chapters, an 
acoustic prediction is carried out over a range dependent 


Simulation to assess how laterally coherent staircases 


impact acoustic propagation. 


B. PROCESS 


The Bellhop model is a beam tracing model for ocean 
acoustic predictions. The theory behind the ray equation 


used in Bellhop can be expressed as: 


= offs) 
d& 1. .ae 
2) See 13 
ds ce or one 
= = ots 
ah 1a 
ds co dz 


where r(s) and z(s) represent ray coordinates in 


cylindrical coordinates, s is the arc length along the ray, 
and c(s) [&s), Cs) ] represents the tangent versor along the 


ray. Initial conditions for r(s), z(s), Cs) and &s) are 


where 0. répresents: Jaunching angle, (f,Z.) 2s: the source 


position, <dand vc, is’ sound speed. at: Source. position. . The 


coordinates are sufficient to obtain the ray travel time: 


ds 
t= [— (15) 
= C(s} 
calculated along the curve [r(s), z(s)] (Rodriguez, 2012). 


To run the desired profiles through the Bellhop model 
an environmental file first needs to be created. Data from 
a C-SALT CTD (Conductivity Temperature Depth) probe from 
station 20 (14°12.70 N 56°17.60 W) that displayed significant 


staircas development was used, as was the MITgcm modeled 


Simulation of laterally coherent well-developed staircases, 


the temperature profile of these examples is in Figure 26. 


The sound speed profile (SSP) was calculated using MATLAB 


code from the seawater toolbox, and is displayed in Figure 


27. The small increase in sound speed at the base of 


individual layers in the staircase region SSP is expected 
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due to the isothermal feature, and increasing depth (and 


hence pressure), of the layers within the staircase. 
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Figure 26. C-SALT Station 20 Temperature profile (left) with 
upper 1000m displayed (right). 
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Figure 27. C-SALT Station 20 SSP (left) with upper 1000 m 
displayed (right). 


The C-SALT SSP was also smoothed to remove the sharp 
changes in sound speed due to the staircase, using a low 
pass moving average filter with a span of 41 points, and 
this profile was used as the example of an = ocean 


environment without staircases, and is in Figure 28 and 29. 
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Figure 28. C-SALT Station 20 SSP (left) and smoothed profile 
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Figure 29. C-SALT Station 20 SSP to depth 700 m (left) and 
smoothed SSP to depth 700 m (right) 
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When building the environment file, only 1000 m of 
data was used, with the bottom type acting as an acoustic 


half-space, with a 1.8 g/cm density and attenuation 


coefficient of 2 dB/m. The attenuation coefficient was 


chosen using Figure 17 of Hamilton (1980), to be just above 


the attenuation value expected for compressional waves 


versus frequency in natural, saturated sediments and 
sedimentary strata, for the frequency range which we are 
focused on (up to 5 kHz). The reasoning behind this was to 
absorb as much sound as possible into the “false” bottom, 
so as to reduce the impact reflected or refracted rays re- 
entering the staircase region from the bottom may have on 


the acoustic propagation. 


The scenarios presented in part b are conducted over 
three source depths: 5 m, 387 m and 450 m. 5 m is a typical 
depth of a fixed hull mounted sonar, and thus provides a 
good representation of a source outside of staircase 
depths. 378 m is the depth at which a rather prominent 
layer interface occurs in the C-SALT profile, and 475 m is 


close to the vertical center of a layer within the 


staircase region. All source depth scenarios were evaluated 


at frequencies of 50 Hz, 2000 Hz, and 5000 Hz to evaluate 
the effect of the deep layers on smaller wavelengths. 
Cz RESULTS 
1. Source Depth 5 m 
a. Ray Trace 


With both the smoothed and stepped profile the 
ray trace (Figure 30) appears the same, regardless of 


frequency. This is largely to be expected as both smoothed 
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and stepped scenarios display the same profile to a depth 


of 250 m, at which point the staircase is introduced. Thus 


a source near the surface experiences surface ducting due 
to the mixed layer, and then a downward refracting ray 
trace due to the decreasing sound speed profile, with a 


shadow zone between the duct and the downward refracted 


rays. The staircase has little impact on this ray trace, as 
the layer depth is not large enough to result in a 
Significant refractive change in the overall downward ray 


trace pattern. 
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Figure 30. Ray Trace profile for a source at 5m ina 
staircase environment. The ray trace is nearly 
identical to that which occurs for a smoothed 

profile. 


b. Transmission Loss 


At a source depth of 5 m, there was no difference 


in the transmission loss between a stepped and smoothed 
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profile. The results were the same as frequency increased. 
Figures 31 and 32 show the transmission loss experienced by 
both profiles at 2000 Hz. This result is contrary to the 
findings by both Wilson (2007) and Wall (2007), which 


stated the source could be outside of the staircase region 


and have a significant effect on propagation. 
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Figure 31. Coherent transmission loss (dB/m) for a source at 
5 m depth and 2000 Hz frequency for a smoothed 
environment. 
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Figure 32. Coherent transmission loss (dB/m) for a source at 
5 m depth and 2000 Hz frequency for a stepped 
environment. 


2. Source Depth 387 m 
a. Ray Trace 


The ray trace for a source at 387 m is 


Significantly different between the smooth and stepped 


profile as seen in Figures 33 and 34. At this depth, the 
source lies at an interface for the stepped profile, but in 
the smoothed profile at this depth the SSP continues to 
decrease. It is evident that when the source lies at the 
same depth as an interface in a region of staircases, there 
is a significant divergence of the rays at this point. For 


low launch angles the rays refract upward or downward, 


leaving a rather significant shadow zone horizontally in 


front of the source. 
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Figure 33. Ray Trace for source at 387 m (depth of an 
interface) in a staircase environment. 
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Figure 34. Ray trace for a source at 387 m in a smoothed 
profile. 
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b. Transmission Loss 


Figures 35 to 40 displays transmission loss for 
both stepped and smooth profiles when the source is 
positioned at an interface, from 50 to 5000 Hz. As was 


evidenced in the ray trace of the stepped profile, a large 


shadow zone exists at this depth at all frequencies tested, 
which is not observed in the smooth profile (although at 


some frequencies there appears to be a small area of 


increased transmission loss above this in the smoothed 


profile). There is also an apparent weak concentration of 


Signal in the area just below this shadow zone, reasons for 


this are unexplained. 
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Figure 35. Coherent TL (dB/m) for stepped profile with 
source in interface at 50 Hz 
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Figure 36. Coherent TL (dB/m) for smooth profile with 50 Hz 
source at 387 m. 
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Figure 37. Coherent TL (dB/m) for stepped profile with 
source in interface at 2000 Hz. 
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Figure 38» Coherent TL (dB/m) for smooth profile with 2000 
Hz source at 387 m. 


Increasing the frequency to 2000 Hz appears to have a 


pronounced effect on the shadow zone in the stepped 


profile. 
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Figure 39. Coherent TL (dB/m) for stepped environment with 
source in interface at 5000 Hz. 
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Figure 40. Coherent TL (dB/m) for 5000 Hz source in smooth 
profile at 387 m. 


There is little obvious difference between 2000 Hz and 
5000 Hz when the source is placed at the interface of a 
staircase 


3% Source Depth 450 m 
a. Ray Trace 


In this scenario the source is positioned at a 
depth that corresponds to the center of a vertically mixed 
layer in the stepped profile. Unlike the previous example, 
in this case the low angle rays converge, causing a very 


obvious duct to occur at the source depth. In the smoothed 


profile, the ducting is not evident to the same extent, 
with no particular features of convergence or divergence 


appearing to stand out. 


51 


Depth (m) 


Sx) 
yo i 
Ret 
ae alte. Ee, 


if 
Nishant 
reread hahah 
a, imatetd 


O 1000 2000 3000 4000 5000 6000. 7000 8000 
Range (rm) 


Ray trace with source at 450 m (center of 


Figure 41. 
in a staircase environment. 


vertically mixed layer) 


KX IK 
OK 
OOK Se 


‘9 ig oe | 
WE 
MG 


= 
= i iB 
a a eAite ty 
8 by a OO, ER ty 
4 fy ; ¥ 
Range (rm) 
Figure 42. Ray trace with source depth 450 m in a smoothed 


profile. 
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b. Transmission Loss 


The ducting that is observed in the ray trace 


using the stepped profile for a source at the center of the 


layer is also observed in all the TL plots. In the smoothed 


profile, all the plots are void of any ducting 


characteristics resultant from staircases. This 


demonstrates that if the source is placed within a layer of 


the staircase, that much higher horizontal acoustic ranges 


may be achieved than if there were no staircase present. 


Once again, examples are shown with lower and higher source 
frequencies. It is evident that an increase in frequency 
from 50 Hz to 2000 Hz results in higher propagation ranges 


within the duct, as a result of increased leakage at lower 


frequencies, as opposed to trapping of smaller wavelengths. 


However, there is not a significant difference between th 


ducted propagation of 2000 and 5000 Hz, and it is proposed 


that this is a result of the increased absorption loss at 


higher frequencies. 
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Figure 43. Coherent TL (dB/m) for stepped profile with 50 Hz 
source in layer. 
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Figure 44. Coherent TL (dB/m) for 50Hz source at 450m ina 
smoothed environment. 
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Figure 45. Coherent TL (dB/m) for stepped profile with 2000 
Hz source in layer. 
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Figure 46. Coherent TL (dB/m) for 2000 Hz source at 450m in 
smoothed profile. 
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Figure 47. Coherent TL (dB/m) for 5000 Hz source in layer of 
stepped profile. 
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Figure 48. Coherent TL (dB/m) for 5000 Hz source at 450 m in 
smoothed profile. 


56 


4. Range-Dependent Source in Layer 
a. Ray Trace 


Following on from the investigation of a source 


placed in a staircase layer as opposed to a staircase 


interface, we investigate a range-dependent xample of 
acoustic propagation in staircases. Figure 24 displayed a 


series of temperature profiles from the numerical model 


spaced every 10 km with lateral coherence over a range of 


300 km. Nine adjacent profiles from this model were 


converted to SSP using the MATLAB acoustic toolbox and then 
compiled into a range-dependent input file compatible with 
Bellhop. Figure 50 shows the individual sound speed 
profiles, and Figure 51 displays an interpolated sound 
speed over the entire range showing the subtle shift in 
sound speed over range and depth as a result of small 


vertical shifts in the laterally coherent layers. The 


center of a 35 m vertical layer at a range of O m within 
the staircase region was identified at a depth of 785 Mm, 
and a ray trace was conducted with the source at this 


depth. Ducting is evident at this depth as can be seen by 


the red ray path in Figure 49. An interesting feature is a 


shift in the ducted rays at a range of approximately 60 km 


from the source. It is hypothesized that this is a result 


of a subtle shift in the depth of interfaces of connecting 


layers as range is increased away from the source. It is 
apparent from the ray trace that the duct shifts up or down 
to remain within the layer, as the layer thickness changes 


(interfaces shift up/down), over extended ranges. 
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Figure 51. Sound speed (m/s) as a function of depth and 
range due to laterally coherent staircases. 


b. Transmission Loss 


Figure 52 to 55 displays transmission loss for 
the range-dependent case when the source is placed at 785 m 
within a prominent layer, over frequencies of 50 to 5000 
Hz. A notable difference in analysis to earlier TL plots is 
that the TL is now being observed over a range of 100 km. 


Ducting is evident at the source depth in all cases with 


the exception of 50 Hz. In the 50 Hz scenario there is 


increased energy at this depth at short range, as a result 


of leakage. However, the larger wavelength means sound is 


not trapped in this layer. The best example of ducting 
occurs at 500 Hz, where there is very obvious trapping of 


the rays through the entire range, with significant energy 
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still remaining at 80 km from the source. The ducting is 


also obvious at 2000 Hz, but is largely diminished by 


5000 Hz. It can be surmised that as the frequency of the 
source is increased the absorption at the higher 
frequencies quickly becomes a more dominant process. 
Therefore, in a range-dependent situation a low to medium 
frequency source of approximately 1000 Hz appears most 
effective in propagation of sound within a duct. It is 
proposed that this frequency is highly dependent on the 
thickness of the vertical layers also, as small layers are 


likely to only trap higher frequency rays. 
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Figure 52. Coherent TL (dB/m) for 50 Hz source within layer. 
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Figure 53. Coherent TL (dB/m) for 500 Hz within layer. 
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Figure 54. Coherent TL (dB/m) for 2000 Hz within layer. 
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Figure 55. Coherent TL (dB/m) for 5000 Hz within a layer. 


5. Range-Dependent Source in Interface 
a. Ray Trace 


Using the same range-dependent profiles from part 
4, a source was placed at a depth of 750 m, which coincides 
with an interface. In the resultant ray trace from this 
there is no apparent ducting of the rays. Similarly to what 


was seen in the range-independent scenario previously 


discussed, there appears to be a divergence of the rays in 
a horizontal direction from the sound source as highlighted 
in the circle. Close to the sound source this may generate 
a shadow zone, however in the range dependent case this 
area is quickly flooded with energy from refracted and 


reflected rays. 
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Figure 56. Ray trace for source at staircase interface for 
range-dependent profiles. 


Transmission loss was again calculated from 50 to 


some energy still manages to get trapped within the lower 


a large 


2000 Hz frequencies. This ducting effect is much less 


pronounced when the source is at the interface as opposed 
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V. CONCLUSIONS 


Thermohaline staircases are a significant feature of 
ocean mixing, and this study has aimed to explore staircase 


evolution in three dimensions as a result of salt 


fingering, through numerical simulation. 


A large task was undertaken in attempting to simulate 


three dimensional staircases over an entire ocean basin. To 


the authors’ best knowledge this is the first basin-scale 


Simulation that exhibits spontaneous layering. A promising 


result has been the modeling of MJlaterally coherent 
staircases over several hundreds of kilometers in the 
ocean. This provides an interesting area of further 
research, particularly with the alignment of layer 
interfaces with background isotherms. While this was not 


achieved in this thesis, it is proposed that Wall’s 


hypothesis of a correlated layer slope with thermohaline 


will hold true, and with more computing hours and time, it 
is believed that a high resolution grid would make this an 


achievable task for researchers in the future. 


This study aimed to parameterize turbulent diffusivity 


through the process of trial and elimination in order to 


closely simulate the thermohaline staircases seen in the C- 


SALT observation area. Through this process, the turbulent 


diffusivity was indirectly estimated at 2.5 x 10° m’/s. 


A comparison between the gamma-effect variation in the 
flux ratio with the density ratio confirmed a very distinct 
link between the flux ratio, the density ratio, and 
Staircase formation. This lends further support to the flux 


ratio model of Radko and Smith (2012). 
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The numerical simulation of merging events indicates 


the prevalence of instability that leads to the ‘*B type 


merger’ as discussed in Radko (2007), and this simulation 


clearly supports the theory that development of 
thermohaline staircases by salt fingers is a direct result 
of strong interfaces growing at the expense of weaker ones 


which ultimately erode. 


The acoustic propagation modeling through the 


staircase of both the C-SALT profiles and the numerical 


model indicate that thermohaline staircases have a very 
pronounced effect on sound propagation, which is 


accentuated with position placement of the sound source 


within the staircase. The ducting that occurs when a source 
is placed centrally within a layer of a staircase is a 
Significant feature, and the increased range of the duct 
with lateral coherence of staircases is an extremely 
important facet to be considered in acoustic modeling. 
Understanding where staircase regions of lateral coherence 
occur, coupled with knowledge of interface and layer depths 
within this region can be of critical importance for 
tactical employment of submarines, gliders, variable depth 


sonars, torpedo homing, and other underwater vehicles 


within this region. This information, combined with an 


awareness of the impact of differing source frequencies on 


the ducting, can result in a significant exploitation of 


the environment for tactical purposes. 
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